Effective field theory: A complete relativistic nuclear model 
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We analyzed the results for finite nuclei and infinite nuclear and neutron matter using the standard 
a — LU model and with the efi'ective field theory. For the first time, we have shown here quantitatively 
that the inclusion of self-interaction of the vector mesons and the cross-interaction of all the mesons 
taken in the theory explain naturally the experimentally observed softness of equation of state 
without loosing the advantages of standard o — u model for finite nuclei. Recent experimental 
observations support our findings and allow us to conclude that without self- and cross-interactions 
the relativistic mean field theory is incomplete. 
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In the quest for an unified model describing both finite 
nuclei and nuclear matter, the quantum hadrodynam- 
ics (QHD) has been a successful tool for the past few 
decades. QHD is the field theory of the nuclear many- 
body problem using hadron degrees of freedom. Models 
based on relativistic QHD takes care ah initio of many 
natural phenomena which are practically absent or have 
to be included in an ad hoc manner in the non-relativistic 
formalism. One of the first successful models based on 
QHD with relativistic mean field (RMF) was constructed 
by Walecka 1] with vector and scalar meson fields. Later 
on, to get reasonable incompressibility and to get good 
results for finite nuclei, cubic and quartic nonlinearities of 
the a meson (standard nonlinear a—oj model) were added 
0. These models were proposed to be renormalizable 
and that constraint limited the scalar interactions to a 
quartic polynomial and disallowed the scalar- vector cross 
interactions and vector-vector self-interactions. However, 
the coupling constants are not assigned with their bare 
(experimental) values but with some effective value to 
have proper results for finite nuclei. Hence the renormal- 
izability of the Lagrangian gets compromised by the use 
of effective coupling constants. 

Inspired by effective field theory (EFT), Furnstahl, 
Serot and Tang abandoned the idea of renormaliz- 
ability and extended the RMF theory by allowing other 
nonlinear scalar-vector and vector-vector interactions in 
addition to tensor couplings H H H H, 0. The EFT 
contains all the non-renormalizable couplings consistent 
with the underlying symmetries of QCD. The effective 
Lagrangian is obtained by employing suitable expansion 
scheme to truncate the infinite number of terms. In such 
scheme the ratios $/Af , W/M, |V$|/M2 and jyVFl/Af ^ 
are the useful expansion parameters |3,y.l5lla,l2| where $ 
and W are scalar and vector meson fields respectively and 
M is the nucleon mass. With the help of the concept of 
naturalness (i.e., all coupHng constants are of the order of 
unity when written in appropriate dimensionless form), 
it is then possible to compute the contributions of the 
different terms in the expansion and to truncate__the _ef- 
fective Lagrangian at a given level of accuracy f" 



None of the couplings should be arbitrarily dropped out 
to the given order without a symmetry argument. Ref- 
erences 0, 0, IS have shown that it suffices to go to 
fourth order in the expansion. At this level one recovers 
the standard nonlinear a — lo model plus a few additional 
couplings, with thirteen free parameters in all. These pa- 
rameters have been fitted (parameter sets Gl and G2) to 
reproduce some observables of magic nuclei 3]. The fits 
display naturalness, and the results are not dominated 
by the last terms retained. This evidence confirms the 
utility of the EFT concepts and justifies the truncation 
of the effective Lagrangian at the first lower orders. 

Recent applications of the models based on EFT in- 
clude studies of pion- nucleus scattering [Ifl] and of the nu- 
clear spin-orbit force j as well as calculations of asym- 
metric nuclear matter at finite temperature with the Gl 
and G2 sets ■ In a previous work Q we have analyzed 
the impact of each one of the new couplings introduced in 
the EFT models on the nuclear matter saturation proper- 
ties and on the nuclear surface properties. In Ref. j9| we 
have looked for constraints on the new parameters by de- 
manding consistency with Dirac-Brueckner-Hartree-Fock 
(DBHF) [H calculations and the properties of finite nu- 
clei. Using EFT we successfully explained the properties 
of the drip- line nuclei as well as the symmetric and asym- 
metric infinite nuclear matter including the neutron star 
and compared with other theoretical calculations p^ . 
Very recently the flow of matter in heavy-ion colli- 
sions is analyzed to determine the pressures attained at 
densities ranging from two to five times the saturation 
density of nuclear matter. This experimental determi- 
nation of the equation of state (EOS) of dense matter 
motivated us to study the applicability of various rela- 
tivistic models at extreme conditions. In this letter we 
present the observations based on the results of our EFT 
calculations, the standard nonlinear a — u model and 
some of the recent experiments. 

The description of EFT and the field equations for 
nuclear matter and finite nuclei can be found in Refs. 

The field equations were derived 0] from an energy 
density functional containing Dirac baryons and classical 
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scalar and vector mesons. According to Refs. 0, 0| the energy density for finite nuclei can be written as 
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where the index a runs over all occupied states (fair) of 
the positive energy spectrum, $ = gs(j)o{r), W = 5„Vb(r), 
R = gpbo(r), A = eAo{r). gs, gv, gp and e are the cou- 
pling constants corresponding to the fields (j)o{r), Vb(r), 
60 (r) and Ao{r) respectively and K3, K4, r]i, r]2, Co, fv 
and fp are non-linear coupling constants. 

The terms with g^, A, Ps and take care of effects re- 
lated with the electromagnetic structure of the pion and 
the nucleon (see Ref. Q). Specifically, the constant g^ 
concerns the coupling of the photon to the pions and the 
nucleons through the exchange of neutral vector mesons. 
The experimental value is 5^/47r = 2.0. The constant 
A is needed to reproduce the magnetic moments of the 
nucleons. It is defined by 
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with Ap = 1.793 and A„ — —1.913 the anomalous mag- 
netic moments of the proton and the neutron, respec- 
tively. The terms with (3s and contribute to the charge 
radii of the nucleon Q . 

Variation of the energy density with respect to (^J, 
and the meson fields gives the Dirac equation fulfilled 
by the nucleons and the meson field equations 0. The 



Dirac equation corresponding to the energy density (fT|l 
and the mean field equations for $, W, R and A can be 
found in Ref. J,, J^J . The meson fields can also be inter- 
preted as Kohn-Sham potentials 0| in the relativistic 
case ^3 aiid ™ this sense they include effects beyond 
the Hartree approach like three-body and many-body in- 
teractions through the nonlinear couplings 0, Q . 

For infinite nuclear matter all of the gradients of the 
fields in the energy density and field equations vanish. 
Due to the fact that the solution of symmetric and asym- 
metric nuclear matter in mean field depends on the ratios 
9s /'^s ^^'^ 9v I ' have seven unknown parame- 
ters. By imposing the values of the saturation density, to- 
tal energy, incompressibility modulus and effective mass, 
we still have three free parameters (the value of g^pjvn^p is 
fixed from the bulk symmetry energy coefficient J). 

The baryon, scalar, isovector, proton and tensor den- 
sities are same as for finite nuclei in the nuclear mat- 
ter limit, i.e., p 
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(Pk = g^fc^ and ^ 
and so on (here kf is the Fermi 



momentum and k is the momentum at any density) . The 
expressions for pressure and energy density are 
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FIG. 1: Scalar (Us) and vector (Uv) potentials against the 
ratio of nuclear matter density with its saturation value (po 
using the parametrizations NL3 0, TMl Gl, and G2 ' 
and in a DBHF calculation 1131. 



Here "f—2 for pure neutron matter and j—A for symmet- 
ric nuclear matter. The asymmetry of nuclear matter is 
defined by the parameter a. For the symmetric matter, 
a=0 and for the neutron matter, a=l. 

While examining the effect of nonlinear coupling in nu- 
clear matter |^ using various nuclear force parameters, 
the NL3 parameter set (considered as a representa- 
tive of the standard nonlinear a — Lu parametrization) has 
been found to fail in following the DBHF results even 
at slightly higher densities. It is well known that the 
DBHF theory in relativistic framework explains well the 
nuclear matter at higher densities 2po) [20I El|. In 
contrast to the standard a — oj model, the EFT calcula- 
tions at high density regimes yield results in accordance 
with DBHF. This scenario is depicted in Fig. 1, where we 
present the results of calculations with TMl parameter 
set [i^ also. In the calculation with TMl, only a quar- 
tic vector self-interaction term is included apart from the 
terms in NL3. This inclusion was done arbitrarily with- 
out considering the underlying QCD symmetries or the 
naturalness. However, with this self-interaction term, the 
TMl give better results at higher densities. This demon- 
strates the importance of self-interactions at higher den- 
sities and exposes the inadequacy of the standard non- 
linear a — uj model. This argument is further supported 
by Fig. 2 in which the variation of binding energy per 
particle (E/A) is plotted as a function of p/po- The NL3 
parameter set gives a much too stiff EOS whereas the 
other parameter sets give a softer EOS which is consis- 
tent with the observed neutron star masses jB'l and radii 
(See Table I) and measurements of kaon production in 
heavy- ion collisions |2^ . 

The recent experimental observations [isl E3 | rule out 
any strongly repulsive nuclear EOS and has confirmed 
the predictions made above. The zero-temperature EOS 
for symmetric nuclear matter derived experimentally is 
shown in Fig. 3 along with the results obtained from 
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FIG. 2: EOS for the same cases as in Fig. 1. 



TABLE I: Upper panel: The surface energy coefficient Eg 
and surface thickness t (in fm). Middle panel: Energy per 
nucleon E/A (in MeV), charge radius r-ch (in fm) and spin- 
orbit splittings Ai5so (in MeV) of the least-bound nucleons. 
Lower panel: the neutron star radius R (in km) and the mass 
ratio M/Mq. 
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Es 




18.51 


18.36 


18.06 


17.8 


16.5-21.0 




t 




1.91 


1.99 


1.98 


2.08 


2.2-2.5 




E/A 




-8.15 


-8.08 


-7.97 


-7.97 


-7.98 




rch 




2.66 


2.73 


2.72 


2.72 


2.73 
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(n,lp) 


5.6 


6.4 
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5.9 


6.2 






(P,lp) 


5.6 


6.3 


5.9 


5.9 


6.3 


«Ca 
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-8.65 


-8.64 


-8.67 


-8.68 


-8.67 
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6.2 
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-8.71 
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-8.71 




Teh 




4.27 
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4.28 


4.28 


4.26 
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1.4 


1.6 


1.8 


1.8 


0.5 
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-7.87 


-7.87 


-7.86 


-7.87 
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5.52 


5.5 


5.5 


5.50 
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(n,3p) 


0.7 


0.8 


0.9 


0.9 


0.9 






(p,2d) 


1.4 


1.6 


1.8 


1.8 


1.3 




R 




13.91 


15.47 


13.93 


10.04 


10.0-12.0 




M/Mq 


2.8 


2.78 


2.16 


2.04 


1.5-2.5 



different calculations. In Fig. 3 we can see that the 
calculations based on NL3 deviate drastically from the 
experimental observation and the EFT calculations with 
Gl upto some extent and G2 to an excellent extent, agree 
with the experiment. A similar situation prevails in the 
EOS of neutron matter and can be seen in Fig. 4. From 
the figures it is very clear that NL3 calculations are not 
suitable for nuclear matter as the formalism is incomplete 
without the self- and cross-interaction terms. On the 
other hand the EFT calculations with Gl and G2 param- 
eter sets explain the situation in nuclear matter naturally 
without any forced changes in parameters or in formal- 
ism and the fit with the experiment is outstanding. The 
experimental data are explained reasonably good only by 
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FIG. 3: Zero temperature EOS for symmetric nuclear matter. 
The shaded region corresponds to experimental data fisl . The 
EOS from calculations using NL3, Gl and G2 parameter sets 
are represented by dashed, dotted and solid lines respectively. 
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FIG. 4: Zero temperature EOS for neutron matter. The up- 
per and lower shaded regions correspond to experimental data 
|15| with strong and weak density dependences, respectively. 
Other cases are same as in Fig. 3. 

the calculations of Akmal el al 23] which employs the 
Argonne Uig interaction. Such an interaction is not ap- 
phed successfully in the case of finite nuclei. The EFT 
calculations are proved to give very good results for fi- 
nite nuclear properties as well as for infinite nuclear 
matter including neutron star . In Table I we present 
some of the sample results of EFT calculations. More 
results for finite nuclei and nuclear matter can be found 
in Ref. Ti"] . In choosing between the parameter sets 
Gl and G2 for further calculations we prefer G2. It is 
worth to note that G2 presents positive values of cou- 
pling constant (^4), as opposed to Gl and to many of 
the most successful RMF parametrizations, such as NL3. 
Actually the negative value of K4 is not acceptable be- 
cause the energy spectrum then has no lower bound |24| . 
However such negative value is necessary in the standard 
a — Lo model to get the results closer to the experimental 
values. On the other hand to have positive value for K4 it 



is not necessary to make two parameter sets as was done 
in Ref. 21]. 

In conclusion, with the new experimental values for 
EOS, the predictions of EFT are proved to be true. With 
the inclusion of self- and cross-interactions and without 
forcing any change of parameters or modifying the for- 
malism the EFT calculations with G2 parameter set ex- 
plain finite nuclei and infinite nuclear matter in a uni- 
fied way with commendable level of accuracy in both 
the cases. Any Lagrangian without all types of self- and 
cross-interactions is incomplete and at present EFT can 
be considered as a complete unified theory for finite nu- 
clei as well as for infinite nuclear matter. 
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